Highly nonlinear dynamics in a slowly sedimenting colloidal gel 
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We use a combination of original light scattering techniques and particles with unique optical 
properties to investigate the behavior of suspensions of attractive colloids under gravitational stress, 
following over time the concentration profile, the velocity profile, and the microscopic dynamics. 
During the compression regime, the sedimentation velocity grows nearly linearly with height, im- 
plying that the gel settling may be fully described by a (time-dependent) strain rate. We find that 
the microscopic dynamics exhibit remarkable scaling properties when time is normalized by strain 
rate, showing that the gel microscopic restructuring is dominated by its macroscopic deformation. 

PACS numbers: 47.57.ef, 64.70.pv, 82.70.Dd 



Gels and attractive glasses resulting from the aggre- 
gation of colloidal particles are the subject of exten- 
sive studies because their physical behavior often results 
from a complex interplay between equilibrium thermody- 
namics and nonequibrium dynamic processes [H-Q, and 
because they are relevant for understanding network- 
forming biological systems Q and for industrial applica- 
tions. Although they exhibit solid-like mechanical prop- 
erties, colloidal gels are easily disrupted by small pertur- 
bations, such as gravitational forces. While a large body 
of macroscopic observations of gels under gravitational 
stress exists [5-13|, very little is known on the micro- 
scopic processes at play during sedimentation, thus lim- 
iting our ability to understand and predict the behavior 
of sedimenting gels. 

Here, we use a novel light scattering method to gain ac- 
cess to the dynamics of a slowly settling colloidal system 
from the macroscopic deformation of the sample down to 
the relaxational behavior at the particle scale. We find 
that the very slow macroscopic deformation occurs via 
irreversible plastic events at the microscopic scale. Re- 
markably, the gel behavior at all scales is controlled by a 
single parameter, the time-dependent compression rate, 
in striking a,nalogy with recent observations on deformed 
polymer and colloidal glasses. 

We study gels formed by attractive colloidal hard 
spheres with radius i? = 82 ± 3 nm, suspended in an 
aqueous solvent at an initial volume fraction i^o = 0.123 
(more details can be found in lB-18|). Gelation is in- 
duced by attractive depletion forces obtained by adding 
micelles of a nonionic surfactant. The interaction be- 
tween colloids is well described by the Asakura-Oosawa 
depletion potential [lij . with a range r « 3 nm. The 
potential can be mapped on the Adhesive Hard Sphere 
model, with a stickiness parameter r ~ 0.01 [ig, [l7j . 
The density mismatch between the particles and the sol- 
vent is Ap = 1.12 g/em'^. The particles have an intrinsic 
optical anisotropy; accordingly, they scatter light with 
polarization both parallel and perpendicular ("depolar- 



ized" ) to that of the incident radiation. The depolarized 
scattered intensity is an accurate probe of the local par- 
ticle concentration [l6l |. 

To probe the sedimentation process in great detail, we 
use a custom-designed light scattering apparatus |19| . 
sketched in Fig. SMI [18j. A low magnification im- 
age of the sample illuminated by a vertical laser sheet 
of thickness 200 /im is formed onto a CCD sensor, us- 
ing depolarized light scattered aX = 90°. By averaging 
the CCD signal over horizontal rows of pixels, we are 
able to measure accurately the time, t, and height, z, 
dependence of the volume fraction, ip{z,t), with a spa- 
tial resolution of about 100 /.tm. We obtain the evolution 
of the local sedimentation velocity profiles, v{z,t), u sing 
an Image Correlation Velocimetry (ICV) algorithm [20[ 
that we apply to rectangular Regions Of Interest (ROIs) 
of the CCD images corresponding to the width of the 
sample and a height Az ~ 0.5 mm. To quantify the mi- 
croscopic dynamics, we measure space and time-resolved 
intensity correlation functions (ICFslfloj. 
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where < • • • >z is an average over all pixels belonging 
to a ROI at height z. Using a mixed spatio-temporal 
correlation method to be described elsewhere, we cor- 
rect the ICFs for the drift due to sedimentation, so that 
(72 — 1 measures the average microscopic particle motion 
between times t and t + t m the co-sedimenting frame 
where v = 0. Because we detect the depolarized scat- 
tered light, 92 — 1 is proportional to the square of the 
self part of the dynamic structure factor, fs{z, q, t, r) = 
exp[—iq-{rj{t+T)—rj {t))], where the sum is over 
the N particles in a ROI at height z and q is the scat- 
tering vector. A full decay of the ICF indicates particle 
rearrangements over a lengthscale = 0.66i?. Thus, 
we are able to measure simultaneously and with both 
temporal and spatial resolution the local volume frac- 
tion, sedimentation velocity and microscopic relaxation 
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FIG. 1: (color online) (a) Temporal evolution of the normal- 
ized height of the gel, h/ho, for cells of section 3x3 mm^ (open 
circles) and 5x5 mm'^ (blue dots). The arrow shows the ap- 
proximate time when the gel is formed. The solid line is the 
prediction of the poroelastic model, Eq. ([l}. Inset: reduced 
height vs. t in a semilog plot. The dotted line is an expo- 
nential fit to the initial decay of the reduced height, the line 
is the fit of the poroelastic model, (b) Concentration profiles 
for various times, as shown by the labels. The lines are fits of 
the poroelastic model, (c) o"(</') obtained experimentally by 
integrating the asymptotic concentration profile {t = 244 h). 
The line is the power law cr[Pa] = 197^9^ (d) Double loga- 
rithmic plot of the velocity profiles v{z) (symbols, labeled by 
time), together with predictions from the poroelastic model 
(solid lines). 



dynamics, in contrast to previous works wore only the to- 
tal gel height, h{t) [i,0,lM{ll, or at most (piz,t) [i,[2lj 
could be measured. 

Figure[T]illustrates the time dependence of h{t), (p{z,t), 
and v{z,t), where z = is the ceU bottom. Up to t « 6 h, 
the sample is partitioned in a supernatant with « at 
the top, an intermediate region with tp = ipo, and a denser 
cake that builds up from the bottom. Rapid fluctuations 
of the scattered intensity indicate that the intermediate 
region is formed by disconnected falling clusters, while 
fluctuations are much slower in the bottom part, reveal- 
ing that a gel has formed. After « 6 h, a relatively 
loose gel {(f < 0.25, Fig. [T}d) occupies the entire colloid- 
rich phase. We follow the slow compaction of the gel 
over more than 10 days, during which its height relaxes 
asymptotically towards a plateau ft-oo ~ 7.3 mm. While 
the time evolution of h(t) ma y b e fitted by an exponen- 
tial decay at short times jlll - ll3l |. deviations are clearly 
observed over days (inset of Fig. [1^), suggesting that sed- 
imentation is controlled by a distribution of relaxation 
times, well reproduced [2^ by the poroelastic model dis- 



cussed below. The total strain of the gel is about 50%, 
such that the gel sediments at an extremely small average 
strain rate of order e « 10~^s~^. As shown in Fig. [TId, 
this very slow macroscopic deformation corresponds to a 
nontrivial inhomogeneous evolution of the local volume 
fraction ip(z,t): ip{z = h{t),t) ?» cpo at all times, while at 
the bottom tp increases up to ip(z = 0) ss 0.5. 

We now turn to the behavior of the local sedimentation 
velocity, a quantity not accessible in previous work. Sur- 
prisingly, the large volume fraction inhomogeneity has no 
effect on v{z, t), which grows linearly with z at all times 
except for the uppermost layer, see Fig. [T]i. This behav- 
ior suggests to characterize the sedimentation with a sin- 
gle time-dependent parameter, the strain rate e = dv/dz, 
a remarkable, but highly nontrivial simplification. 

Our detailed set of measurements allows us to perform 
a rigorous quantitative test of the modeling of gels as a 
"poroelastic" medium pionee red by B uscall and White Q 

The modeling 
Our velocime- 



and used in several studies [il-13, 21 
simplifies in the absence of wall friction 
try analysis indicates that there is no horizontal compo- 
nent of the displacement resulting from the compression 
of the gel, i.e. that the (effective) Poisson ratio of the gel 
is negligible, a remarkable property that was discussed 
before |ll| but could not be measured experimentally. 
Accordingly, no significant fraction of the gravitational 
stress is redirected on the walls, making wall friction 
vanish, as suggested also by the observation that chang- 
ing the cell section from 3x3 mm^ to 5 x 5 mm^ does 
not change the evolution of /i, ip, v nor that of the mi- 
croscopic dynamics discussed below. Neglecting inertial 
terms, by balancing the gravitational stress on a gel slice 
by the sum of the viscous stress, due to the flow of the 
suspending fluid through the network, and the restoring 
stress resul ting from the gel response to a deformation, 
one gets [1, [23| : 
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where K(ip) = ipda/df is the effective compressional 
modulus in response to an applied stress a, and k{(p) is 
the gel permeability, a constant on the order of the typ- 
ical pore cross section that relates the solvent flow rate 
to the pressure drop in Darcy's law; rj and g are the sol- 
vent viscosity and the gravitational acceleration, respec- 
tively. The boundary conditions are ip{z = h,t) ~ ipo 
and zero mass flux at the cell bottom. The initial con- 
ditions are obtained from the concentration and velocity 
profiles measured experimentally at t = 16 h. 

The solution of Eq. ([T]) depends crucially on the ma- 
terial properties via K{ip) and K{ip). The compressional 
modulus may be estimated from the concentration profile 
at large t, when v k, Q , Il7l |. For t — ^ oo, the viscous 
stress vanishes and the restoring stress tj[(p(z)] at a given 
height balances the total weight per unit area of the par- 
ticles lying above. Thus, a numerical integration of the 
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equilibrium profile along z yields u{}p) [13|,|l7[. As shown 
in Fig. [TJ;, this yield stress is well fitted by cr = ayj'^^^ '^, 
where the exponent is close to that measured for the de- 
pendence of the linear elastic shear modulus of colloidal 
gels with short-ranged interactions [24 1. 

Various forms have been proposed in the past for the 
permeability ^^{^p)^ but no stringent tests on the validity 
of the competing expressions could be performed, owing 
to the lack of detailed measurements of the sedimenta- 
tion kinetics. Expressions assuming a fractal morphology 
of the gel fail to reproduce our data, presumably be- 
cause the volume fraction of our gel is too large for a 
fully fractal structure to develop. Instead, we find that 



(1-y)" 



as in Ref. [25|, reproduces very well 



the behavior of h{t), see Fig. [T^. The fitting parameters 
are kq = 6.78 x 10^^ A'm^, corresponding to a typical 
pore size in the range OAR — 5.9R (depending on cp), 
and m = 7, slightly higher that m = 5.5 — 6.5 typi- 
cally reported for hard spheres [H, [l^. With kq and 
m fixed by the fit to h{t), Eq. ([T]) is solved to obtain 
the temporal evolution of the full concentration and ve- 
locity profiles. Remarkably, we find that the poroelas- 
tic model captures very well -with no further adjustable 
parameters- both ip{z, t) and v{z, t), as shown in Figs.[TjD 
and Id, demonstrating its efficient modeling of the sed- 
imentation kinetics of our gel at both macroscopic and 
mesoscopic scales, and suggesting that the sedimentation 
occurs slow enough for the effective compressional mod- 
ulus to approach its asymptotic limit at e — !• 0. 

Although very successful, the poroelastic model pro- 
vides no insight on the microscopic dynamics and its 
relation with the macroscopic deformation, a key issue 
unexplored so far. Given that K{(p) appears to be unaf- 
fected by the settling kinetics, one wonders whether the 
microscopic dynamics of the gel is also independent of 
sedimentation. Indeed, it has been argued that gravita- 
tional compaction provides a convenient way to slowly 
change ip in order to probe the equilibrium dynamics 
of attractive cqUoidal systems as a function of particle 
concentration [2(||, implicitly assuming that microscopic 
dynamics and sedimentation decouple. 

To investigate the link between microscopic rearrange- 
ments and the macroscopic gel compaction, we examine 
the intensity correlation functions g2iz, t, t) — 1. Figure[2] 
shows (72 — 1 for several heights and times during the gel 
sedimentation. At all i, a full decay of (72 — 1 is observed, 
indicating that the microstructure is continuously chang- 
ing. At large t, an intermediate plateau becomes increas- 
ingly visible and 52 — 1 exhibits a two-step decay. We re- 
port in Fig.[2}3 the decay time of the final relaxation of the 
ICE, Ti/e, defined by 52(^/6) - 1 = e~^[g2{T = 1 s) - 1]. 
At any given time, Ti/e is nearly independent of z. Thus, 
the microscopic dynamics is fairly homogeneous across 
the sample, in striking contrast with the highly inhomo- 
geneous concentration profiles reported in Eig.[T|D. While 
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FIG. 2: (color online) Main panel: intensity autocorrelation 
function measured at different sedimentation times as indi- 
cated by the labels, for a cell of section 3x3 mm^ . To avoid 
overcrowding the graph, we plot only data for a representative 
height, 2 = 5 mm. Inset: relaxation time ri/e as a function 
of height 2 for various times (same symbols as in the main 
plot). For all times, Ti/e does not depend on z. 



volume fraction usually controls the spontaneous dynam- 
ics occurring in attractive gels and glasses 0,131, we con- 
clude that ip plays no such direct role in our system. 
Similarly, Ti/^ is independent of froi = ^^/(l ~ ths 
fluid velocity relative to the particles, because Uroi grows 
with z (data not shown) while ri/e is z-independent. This 
rules out the possibility that the dynamics is due to re- 
arrangements induced by the flow of solvent through the 
gel structure. This is further substantiated by the fact 
that the strain rate associated with the gel deformation 
(e < 10~^ s""'^, see Fig. [3]) is at least two orders of mag- 
nitude smaller than that sustained without flow-induced 
damage by gels in typical linear rheology experiments. 
Finally, the decay of 32 — 1 cannot be due to the relative 
motion of particles imposed by the affine vertical com- 
pression of the gel, since the ICE is sensitive only to the 
horizontal component of the displacement, parallel to q. 

Instead, we suggest that structural reorganization is in 
fact a highly nonlinear process induced by the very small, 
but finite, strain rate e of the sedimenting gel. In col- 
loidal [Til, polymeric 14 1 and molecular 27 1 glasses un- 
der stress, the microscopic dynamics is unaffected by the 
fiow when the inverse strain rate is much larger than the 
equilibrium structural relaxation time. Because the lat- 
ter can become very large in amorphous materials, very 
small rates can enhance particle mobility, to the point 
where the macroscopic fiow rate controls the microscopic 
relaxation time. To test this hypothesis, we plot in Fig. [3] 
Ti/o vs. the measured strain rate e. The data are very 
well fitted by a power law, ti/^ ~ ^-o.98±o.02 ^^gj. jjjQj.g 
than two orders of magnitude, establishing a direct re- 
lationship between microscopic structural relaxation and 



macroscopic flow, Ti/c ^ ^ , analogous to that observed 
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FIG. 3: (color online) Main panel: double logarithmic plot of 
Ti/c us. e. Solid line: power law fit yielding Ti/g ~ ^-o.98±o.o2 
The solid symbols are as in Fig. [51 the open symbols are data 
for a cell of section 5x5 mm^ and time t = 18, 25, and 30 h 
for squares, circles and triangles, respectively. Inset: scaling 
of 32 — 1 measured at various heights and times (see labels, 
where height is in mm), when plotted against er. The line is 
the function 0.35 exp(— 75er). 



in strongly sheared glassy materials. As a further test, 
we plot in the inset of Fig. [3] ff2 — 1 against rescaled time 
CT. Although the initial decay of 172 — 1 varies with t, the 
second step of the relaxation of the ICFs measured at dif- 
ferent heights and times collapses reasonably well onto a 
master curve (line). Note that structural rearrangement 
due to plastic events occurs for a cumulated strain eri/o 
of the order of 1%, comparable to typical yield strains 
measured in colloidal systems [lij or molecular glass for- 
mers [2i,[2i|. 

Our experiments provide direct evidence that in a col- 
loidal gel sedimenting over several days, the small ex- 
ternal perturbation induced by gravitational forces has 
dramatic consequences at the particle scale, where struc- 
tural relaxation continuously occurs at a rate imposed 
by the macroscopic deformation of the sample, in strik- 
ing analogy with the enhancement of particle mobility 
observed in deformed repulsive colloidal glasses [l^ and 
polymeric glasses [14| | . Given the intricate and non- linear 
nature of the microscopic dynamics, it is all the more re- 
markable that the simple macroscopic poroelastic model 
widely used in sedimentation and filtration problems can 
successfully describe the sedimentation kinetics of our gel 
to a remarkable accuracy, by specifying only a couple of 
effective mechanical properties of the material. 
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